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ABSTRACT 

We present Spitzer infrared (IR) photometry and spectroscopy of the lensed Lyman break galaxy 
(LBG), MS1512-cB58 a,t z — 2.73. The large (factor ~ 30) magnification allows for the most detailed 
infrared study of an = 3) LBG to date. Broadband photometry with IRAC (3-10 /um), IRS 

(16 /im), and MIPS (24, 70 & 160 /im) was obtained as well as IRS spectroscopy spanning 5.5-35 /im. 
A fit of stellar population models to the optical/near-IR/IRAC photometry gives a young age (~ 9 
Myr), forrming stars at ~ 98 M0 yr"-'^, with a total stellar mass of ~ 10^ Mq formed thus far. The 
existence of an old stellar population with twice the stellar mass can not be ruled out. IR spectral 
energy distribution fits to the 24 and 70 fim photometry, as well as previously obtained submm/mm, 
data give an intrinsic IR luminosity i/^ = 1 — 2 x 10^^ Lq and a star formation rate, SFR ^20 — 40 
M0 yr~^. The UV derived star formation rate (SFR) is ~ 3 — 5 times higher than the SFR determined 
using Lm or Lna because the red UV spectral slope is significantly over predicting the level of dust 
extinction. This suggests that the assumed Calzetti starburst obscuration law may not be valid for 
young LBGs. We detect strong line emission from Polycyclic Aromatic Hydrocarbons (PAHs) at 6.2, 
7.7, and 8.6 fim. The line ratios are consistent with ratios observed in both local and high redshift 
starbursts. Both the PAH and rest-frame 8 /im luminosities predict the total i/fl, based on previously 
measured relations in starbursts. Finally, we do not detect the 3.3 /im PAH feature. This is marginally 
inconsistent with some PAH emission models, but still consistent with PAH ratios measured in many 
local star-forming galaxies. 

Subject headings: galaxies: high-redshift, galaxies: individual (MS1512-cB58), infrared: galaxies 



1. INTRODUCTION 

Much of the global star formation at z > 2 occurs 
in ultraviolet -luminous Lyman Break Galaxies (LBGs, 
ISteidel et al!l [r996: Reddv et al.l l2005l. For most LBGs, 
much of what we know about their star-formation (star 
formation rates, dust reddening) is derived from rest- 
frame UV properties, where considerable degeneracies 
exist (ie. dust reddening and starburst age). In most 
L^jy LB Gs, the majority of the UV p hotons are absorbed 
by dust (|Adelberger fc Steidell 120001 : iReddv et al.ll2008D . 
which then emits the energy in the infrared. Therefore, 
an accurate measurement of the total infrared luminosi- 
ties would give a complete census of the reprocessed UV 
photons and thus, a better determination of the star for- 
mation rates in LBGs. LBGs can be detected in the 
infrared at 24 /zm with the Spitzer Space Telescope. Ob- 
servations in the Spitzer 24 /im bandpass are the most 
sensitive towards detecting dust obscured star formation 
at 2; < 3. However, even the deepest Spitzer surveys 
(/24 > 20 /iJy, or L/_r^3 x 10^^ L© at z = 3), can not 
de tect the majority o f LBGs (see IR luminosity function 
of IReddv et"aIll2008^ . Between 1 < z < 3, the 24 /im 
flux is dominated by the emission features of Polycyclic 
Aromatic Hydrocarbons (PAHs), resulting in uncertain 
bolomctric corrections due to variations in the PAH emis- 
sion and possible silicate absorption at 9.7 fim. As the IR 
spectra of typical LBGs have not been measured, tem- 
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plates of local starbursts have been used to extrapolate 
the rest-frame 5-12 /xm luminosities to total IR luminosi- 
ties. It is not clear, however, that the IR SEDs, and in 
particular the PAH emission, is the same in local and 
high redshift starbursts. 

Much progress has been made in measuring the 
mid-IR spectral properties and broadband SEDs of 

dus ty, ultr aluminous star- f orming gal axies at z > 

1 (lYan et al. 200 51: iMenen dcz-Delmcst re et all 120071 : 
ISaiina et all i2007i : iDesai et al.i i2007i : iPope et al.ll2008f l. 
and even an ext remely luminous and dusty LBG 
(jHuang et al.ll2007l ). However, there have been no mid- 
IR spectra obtained of LBGs with typical luminosities 
(Afisoo ^ —21, i/_R, ^ 10^^ Lq) and dust extinction 
(0.0 < E{B - V) < 0.3). There are a few LBGs 
that are sufficiently magnified through lensing for their 
mid-IR spectral pr operties and far - IR photome t ry to 
be measured (e.g ., lYee et al.l Il996t ISmail et al.l 120071: 
lAllam et al.|[2007f) . MS1512 -cB58 (hereafter , cB58) is the 
first lensed LBG (z ^ 2.73. lYee etaD '1996') to be found 
with a large m agnification 30, Williams & Lewii 
Il996t ISeitz et ar.lll998f ). allowing extensive follow-up at 
various wavelengths. Fits of stellar population models to 
the optical and near-IR photometry (rest-frame UV and 
optical) indicate a very young {tage ~ 10 — 20 Myr), rela- 
tively dusty, E(B — V) ^ -3, starburst (jEllingson et all 
\199&} . iPettini et all (|2000( l use the rest-frame UV con- 
tinuum to determine the dust extinction and derive a 
reddening (and magnification) corrected star formation 
rate of 93 Mq yr~^. A Keck near-IR spectrum of nebu- 
lar emission lines y ields an intrinsic st ar formation rate 



of ~ 20 Mp 



(iTeplitz et al.ll2000l) Sub-millimeter 
van der Werf et al. 1 120011 : IS awic kil l2001h and millimeter 



Baker et al.ll2001l )~photometrv also suggest, though with 
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large uncertainty, a star formation rate significantly lower 
than the UV-derived SFR, calling into question the ap- 
plicability of locally derived relations (reddening laws, 
IR SEDs) to this LBG, and perhaps the population as 
a whole. However, this discrepancy can be reconciled 
if the IR SED is very w arm, as predicted by th e LBG 
dust emission models of lTakeuchi fc Ishiil (|2004D . Here 
we present Spitzer mid- and far-IR photometry, and mid- 
IR spectroscopy of cB58 to more accurately determine its 
IR properties and whether local starburst relations can 
be used for analysis of LBGs. 

We use a ACDM cosmology with = 0.3, flA — 0.7, 
and Hq = 70 km s~^ Mpc~^. All intrinsic luminosities 
and star-formation rates are corrected assu mi ng a lensing 
magnification, /i = 30 (jSeitz et al.l ^1998V Scit z et all 
([TQQSi ) estimate an error on the magnification 30%. 

2. OBSERVATIONS AND DATA REDUCTION 

Both IRAC and MIPS observations were part of a GTO 
program (PID:65, PI: C. Lawrence). cB58 was imaged in 
aU fo ur bands of the I RAC instrument (3.6/4.5/5.8/8.0 
/xm, iFazio et aLll2004[ ) 2004 January 18. In each band, 
the target was imaged at nine dither positions for a to- 
tal of 900s per ba nd. Far-IR photometry with MIPS 
(jRieke et al.l 120041) was taken 2004 March 20. The 24, 
70, and 160 /xm bands were imaged for 2263, 1091, and 
105 seconds, respectively. 

All IRS observations were taken for Spitzer Program 
30832. We obtained 16 /im with the Blue Peak-Up im- 
ager on the Spitzer IRS instrument on 2007 June 25 for a 
total exposure time of 1415 seconds. IRS spectral obser- 
vations were performed on 2007 March 23 (Long-Low) 
and 2007 June 21 (Short-Low). The spectroscopy was 
taken in mapping mode, placing the target in four differ- 
ent positions along the Short-Low slit and six positions 
along the Long-Low slit as recommended bv lTeplitz et al.l 
(|2007l ) for optimal S/N in deep IRS integrations. All IRS 
observations were scheduled immediately after the "sky- 
dark" calibrations to reduce the chance of latent images 
from preceding bright targets. Total integration time was 
6.77 ks for both SLl and SL2, 14.63 ks in LL2, and 34.13 
ks in LLl. 

cB58 is only 5.3" away from the cD galaxy of the 
foreground cluster {z = 0.373) and this galaxy has ap- 
proximately the same flux at 24 /im (250 /iJy) as cB58, 
so extra caution was required to avoid light from the 
cD galaxy falling in the IRS slits. For both Short-Low 
and Long-Low observations, the slit was aligned as close 
to perpendicular with the cB58-cD connector as was al- 
lowed given scheduling constraints (see Figure [J). Since 
the Long-Low slit is 10.7" wide (approximately twice the 
distance between cB58 and the cD galaxy), half of the 
light from the cD galaxy would have fallen into the slit as 
well. Therefore, we offset the IRS Long-Low observations 
2.35" away from the cD galaxy. Given the point-spread 
function, we estimate that this eliminates ^ 87% of the 
cD galaxy light while only decreasing the cB58 flux by 
-6%. 

The IRAC data were reduced using the SSC's MOPEX 
software and drizzled on to a finer pixel scale (0.6" 
pix~^). The cD galaxy is bright in the IRAC images 
and it's profile clearly extends beyond cB 58 (Figure [1]). 
We use the GALFIT package ()Peng et al.i,2002.) to fit a 
de Vaucouleurs profile to the cD galaxy and subtract it 



from the images. We then extract the flux in a large el- 
liptical (8"x6") aperture (cB58 is marginally extended). 
We do not apply aperture corrections because these are 
small (< 3%). At 3.6 and 4.5 fim, where the cD galaxy is 
brightest, GALFIT has slightly oversubtracted the pro- 
file of the cD galaxy in the region of cB58. Therefore 
the background is slightly lower around cB58. We add 
a small constant to the background to correct for this, 
resulting in a ~ 5% increase in the cB58 fluxes at 3.6 & 
4.5 /im. We added a 5% systematic error (to account for 
the uncertain background estimation) in quadrature with 
the statistical errors and the uncertainty in the absolute 
IRAC calibration (5%). 

The IRS 16 /im and MIPS 24 /im data were reduced 
with the MOPEX package. The 70 /im data were re- 
duced using the Germanium Reprocessing Tools (GeRT), 
following the t echniques opt imized for deep photometry 
data given by iFraver et all (|2006) . The IRS 16, and 
MIPS 24 & 70 /im images are shown in Figure [1] The 
PSF at 16 /im is sufficiently small to resolve cB58 from 
the cD galaxy. In both the 24 and 70 /im images, confu- 
sion with nearby sources is a concern. At 24 /xm, cB58 is 
marginally confused with the cD galaxy, 5.3" to the SE 
(the 24 /im FWHM is 6"). We use the IRAC positions as 
priors to simultaneously fit two point sources at these po- 
sitions, giving /24(c-B58) = 240 ±40 /iJy. At 70 /im, the 
confusion is more severe, as there is also a bright 24/70 
micron source ~ 12" N-NW of cB58 (the 70 /im FWHM 
is 18"). Three point sources are fit simultaneously at the 
24 /im positions, giving a /7o(ci358) — 1.7 ± 1.0 mJy. 
The cD and the northern galaxy are both brighter (2.7 
and 4.5 mJy, respectively). The error includes the error 
in the multiple source fit. If we do only a two compo- 
nent fit with the bright NW galaxy and the combined 
cB58-|-cD, we get a fjo{cBb9> + cD) = 4.7 mJy, a con- 
servative upper limit to the cB58 flux. We do not detect 
cB58 at 160 /im and derive a 3(t upp er limit of f iRn < 24 
mJy based on the confusion limit of iDole et al.l ()2004[ ). 

IRS data r educt ion was performed as specified in 
iTeplitz et all (l2007t l. First, we removed latent charge 
by fitting the slope of the increase in backround with 
time and subtracting this background row by row. Sec- 
ond, "rogue" pixels were masked using the IRSCLEAN 
program provided by the SSC. Finally, the observations 
at other map positions were used to determine the sky, 
which was then subtracted. The individual frames were 
co-added to produce 2D spectra at each map position. 
One-dimensional spectra were optimally extracted at 
each map position using the SPICE software provided 
by the SSC. 

As mentioned above, special consideration is required 
to account for contamination of the cB58 spectrum by the 
cluster cD galaxy. Contamination is primarily of concern 
at longer wavelengths (Long-Low spectrum) because of 
the broader PSF, and the wider slit width (10.7"). We 
estimate the contribution of the cD galaxy to the IRS 
Long-Low spectrum in the following manner. First, we 
use the SSC's model of the 24 /im PSF (for a T = 50K 
blackbody spectrum) and determine the fraction of the 
flux falling into the Long-Low slit for several slit posi- 
tions. For a slit offset from the center of the PSF by 
2.35" (the Long-Low position offset from cB58), we de- 
tect ^ 93.6% of the flux that is detected when the source 
is centered within the slit. For a slit offset 7.35" from the 
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PSF center (the total offset from the central cD galaxy), 
we detect only ~ 23.4% of the flux detected when the 
source is centered within the slit. Because the 24 //m 
fluxes of the two sources are known (/24(cB58) = 240 
& f2i{cD) — 250 /iJy), we determine that the resulting 
spectrum will have a contribution of 250 x 0.234 = 59 
/xJy from the cD galaxy and 240 x 0.936 = 222 /itJy from 
cB58 for a total of 281 /xJy. This agrees very well with 
the extracted spectrum flux (uncorrected for contamina- 
tion) /24 = 284 ^Jy, when multiplying by the MIPS 24 
/im transmission curve and integrating over the band- 
pass. Therefore, to correct the spectrum for contam- 
ination from the cD galaxy, we subtract 23.4% of the 
estimated cD galaxy spectrum. Then, to correct for slit 
loss of the cB58 flux we multiply the spectrum by 1.064. 

UicBbS) = 1.064[/^(o6s) - 0.2MU{cD)] (1) 

The cD galaxy spectrum, /^(cD), is assumed to be 
a power law flt to the 16 & 24 /im fluxes because, 
at the redshift of the cD galaxy {z = 0.37), the IRS 
Long-Low 1st order spectrum probes wavelengths (16 < 
Arest < 24 ^m) where strong emission lines or absorp- 
tion features are not typically seen in either star- forming 
galaxies or AGN (lArmus_etk|[20ql iBrandl et"aH l2006t 
iBuchanan et al.ll2006l: iNetzer et al.l[2007f ). The net effect 
is that the continuum is supressed by 10-20% while the 
PAH fluxes remain largely unchanged. 

This exercise only tells us the fractional contamination 
levels at ~ 24 /im, but the contamination is a function 
of wavelength because the PSF width is a function of 
wavelength. The contamination fraction was computed 
over the large 24 /im bandpass (21-27 fj.Jy) so the PSF 
will not be significantly different between the 20-35 /im 
region of interest. Therefore we do not estimate a wave- 
length dependent contamination correction as this is a 
second order correction. 

3. RESULTS 

3.1. Stellar Population Fits 

The Spitzer photometry of cB58 is given in Table [TJ 
The IRAC photometry traces the rest-frame near-IR, al- 
lowing us to more accurately determine the stellar popu- 
lation parameters than with rest-frame UV/optical data 
alone. We combine the o ptical an d near-IR photome - 
try from lYee et all (|1996f ) and EUi ngson etldl (|1996D . 
respectively, with the IRAC photometry to fit the com- 
plete stellar SED from rest-frame 1200 A to 2 Aim. We 
use population synthesis models from iBruzuall (|2007D . 
which account for the thermally pulsing AGB phase of 
stellar evolution. Three different star formation histories 
(constant, single burst, exponential decay with e-folding 
time, T — 100 Myr) are used and the age of the burst is 
allowed t o vary. The du st obscuration la w for starburst 
galaxies (|Calzettilll997l) . a Salpter IMF (|Salpeteilll955f) 
and a metallicity Z = OAZq is assum ed, similar to previ- 
ously deter mined cB58 metallicities ([Pettini et al.ll2000l : 
iTeplitz e~a l. 2000). We find that a very young, con- 
stant star formation model is required to explain the 
lack of a Balmer break between the J and H bands. 
The best-fit SED, shown in Figure [H has an age of 
9.313;} Myr, SFR = 98^3^ Mq yr'^, an above aver- 
age reddening for an LBG, E{B ~ V) ^ 0.35 ± 0.03, 
and a total stellar mass created in the current episode of 



log(M[Mo])= 8.9 4 ± 0.15. As has be en found in other 
studies of LBGs (jShaplev et al.l [20051 ). the inclusion of 
IRAC photometry gives stellar population parameters 
similar to those derived with fits to op tical/near-IR pho- 
tometry alone (|Ellingson et al.|[l996f ). but with smaller 
uncertainties. 

With the addition of the rest-frame near-IR (IRAC) 
data, it is possible to look for a population of older, 
redder stars that would boost the IRAC photometry 
above what is e xpect ed from the young starburst alone. 
iBechtold et all (| 19981 ) use ISO mid-IR (6.7 & 11.5 /im) 
photometry and conclude that no more than 10% of the 
mass is from an old (~ 2 Gyr) population. However, be- 
cause the starburst in this system is ~10'^ years old, the 
rest-frame near-IR is highly luminous due to red super- 
giants. This significantly decreases the near-IR mass-to- 
light ratio relative to models at other ages (eg. 10^ or 10^ 
yrs). Thus, a massive older stellar population can exist 
but negligibly affect the IRAC photometry. We show an 
example in Figure [21 where we include an older (2 Gyr 
instantaneous burst) stellar population with more than 
two times the stellar mass of the younger stars (~ 2 x 10® 
Mq). The best fit reddening of the older stellar popu- 
lation, E{B — V)oid = 0.27, is somewhat less than the 
younger burst component, E{B — V)young — 0.38. The 
rest-frame near-IR flux of the older population is still a 
factor of two lower than that of the less massive, younger 
population. In this situation, the addition of the two 
components has nearly the same SED as the single com- 
ponent flt. We conclude that ~ 10® Mq has been pro- 
duced by the current star formation episode thus far, 
but there can be twice that mass in older stars while still 
fltting the data reasonably well (reduced = 2.5 vs. 
1.7 for single component flt), resu lting in a maximum 
Mstars < 3 X 10® Mq. Howcvcr, iPettini et all (|2002f ) 
show that the lack of N and Fe-peak elements in the ISM 
suggests that most of the stellar mass is less than 300 
Myr old. Therefore, it is likely that the mass from an old 
stellar population is signiflcantly less than the maximum 
amount obtained from our two component flt. (Note: 
The stated masses are the current existing stellar mass 
and do not include the mass returned to the ISM via 
SNe. For the old (2 Gyr) single burst model, the original 
total mass is equal to 1.37 times our stated values. The 
correction to the young, '^lO Myr models i s negligible.) 

Th e velocity dispersions in atomic (H a, [ Teplitz et al.l 
I2000f ) and molecular (CO, iBaker et aL 2004) _gas both 
give virial masses of ~ IQio Mq. [Baker et all (|200i ) 
derive a molecular gas mass of Mgas = 6.6^4 3 x 10® Mq. 
Therefore, our maximum stellar mass plus molecular gas 
mass does not exceed the maximum total mass. 



3.2. Far-Infrared Spectral Energy Distribution 

The MIPS photometry is plotted in Figure [3] We plot 
both the conservative 70 /tm upper limit as well as the 
valu e derived from the three component fl t. The 850 
jiui van der Werf et aI1[200l[) and 1.2 mm Baker et al.l 



I2OOII ) detections are plotted as well. Note that a 3a up- 
)er limit of 3. 9 mJy at 850 pm has also been reported 
Sawickill2001[ ). in disagreement wit h the nearly 5g re- 
orted detection fs5o = 4.2 mJy bv Ivan der Werf et al.l 
2001f ). To estimate the IR lum inosit y, we flt the 105 in- 
frared SEDs of Charv fc Elbai pOOll . hereafter CEOl) to 
the detections at 24, 70, 850, and 1200 /^m. The best-flt 
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model is shown as the sohd curve in Figure [3l The total 
infrared (8-1000 /im) luminosity is Lm = 1.1 ±0.1 x 10" 
Lq. If we assume that the local templates of CEOl ap- 
ply at high-z, and we simply use the template with the 
same L24, we get an Lm — 1.5 x 10^^ L0 (dashed line 
in Figure [3|) . This is actually the same template SED 
as in the complete SED fit, but with a higher normal- 
ization. We also fit the SEDs of iDale fc Heloul (|2002l 
hereafter DH02) and get a best-fit L//? = 2.4 x 10" Lq, 
about two times higher than with the fit to the CEOl 
templates. However, the DH02 fit is significantly worse 
than the best-fit CEOl models (reduced = 4.8 with 
DH02 vs. 1.6 with CEOl). The DH02 models do not vary 
significantly in the predicted submm-to-24 /im flux ratio 
and every model underpredicts the /24 given the 850 /im 
& 1.2 mm fluxes. We therefore adopt the best-flt CEOl 
model and IR-luminosity but estimate a significant un- 
certainty of ~ 0.3 dex. The Lm < 2.4 x 10^^ L© is a 
reasonable upper limit since the best-fit DH02 template 
passes just under our 3cr 70 & 160 /im upper limits. 

Finally we note that, although the CEOl mod- 
els produce better fits than the DH02 models, none 
of these models has the large observed /24//850 or 
/24//1200 ratios observed in cB58. The observed 
/24//85o(ci?58) = 0.057 is larger than the maxi- 
mum template values f24:/ f85o{CE01)max — 0.055 
and f2i/fs5QiDH02)macc = 0.034. The observed 
/24//i20o(ci?58) = 0.23 is larger than the maxi- 
mum template values f2i/ fi20o{C E01)max = 0.13 and 
f24/ fi2oo{DH02)„iax = 0.06. There are starbursts at 
low redshift that have higher /24//850 ratios than cB58. 
M82, for example, would have a ratio /24 //850 = 0.08 at 
z = 2. 73, when using the compiled SED of lGalliano et al.l 
()2007t ). Therefore, the models simply do not display the 
full range of variability in the SEDs of local starbursts. 

3.3. Mid-Infrared Spectrum 

The mid-infrared spectrum of cB58 is plotted in Fig- 
ure U) As a reference, we have also plotted the Long-Low 
first order spectrum before the corrections for contamina- 
tion and slit loss explained in Section [2l The corrections 
make only small differences to the overall spectrum be- 
tween 20 < A < 33 /im. Strong PAH emission is seen at 
6.2, 7.7, and 8.6 /im. The line strengths were measured 
by simultaneously fitti ng Prude profile s with centers and 
widths deflned by Drai ne fc ll ()2007( ). We also flt an 
underlying continuum and allow both the slope and am- 
phtude to vary. The resulting flt is plotted in Figure El 
There is clearly large uncertainty in the amplitude of the 
minor PAHs, but the fluxes of the stronger PAHs (6.2, 7.7 
complex, and 8.6) are well constrained. The PAH fluxes 
are listed in Table [3l The continuum component is very 
low (the best-fit continuum comprises only ~ 13% of the 
total flux measured between 5.7 < Xrest < 9.0/im), and 
is not required for a good flt to the spectrum. However, 
because we do not have data at Xrest > 9/im, we can not 
completely rule out a somewhat larger continuum compo- 
nent combined with silicate absorption at Xrest = 9.7/im. 

We do not detect the PAH emission feature at A = 3.3 
/im, putting a limit on the PAH ratio of Lq ^/ L^,^ > 3.9. 
This ratio depends strongly on the i onization state and 
gram size di stribu tion of the PAHs (|Li fc Drainell2001t 
iDraine fc Lil l2007f ) . Assuming the grain size distribu- 
tion and ionization fraction fit for Milky Way dust, the 



iDraine fc Lil ()2007f ) models predict a luminosity ratio 
-^6.2/^3.3 ^ 3, suggestin g that we sh o uld s ee the line 
at ~ 4cr. We note that IDraine fc Lil ()2007D increased 
the strength of the 3.3 /i m feature by a fact or of 1.5-2.0 
higher than the values in 'Li fc Draind ()2001l ) to better fit 
theoretical calc ulations of the abso rption cross section for 
various PAHs (|Malloci et al.ll200l . 

The best empirical dataset for fl ux ratios b etween the 
6.2 and 3.3 /im featu res is that of llmanishil ((2006) and 
llmanishi et al.l ()2007l ). where 6.2 /tm is measured by 
Spitzer and 3.3 /im from the ground. The sample con- 
tains nine ULIRGs that have no strong evidence (radio, 
optical, mid-IR) for an AGN and are compact enough 
(ie. no merging components) such that varying slit loss 
is a smaller issue. These nine ULIRGs, have an average 
< ^6.2/^3.3 >= 5-6, with 7/9 having a ^6.2/^3.3 larger 
than our 3ct limit. Thus, if the PAH ratios in LBGs are 
similar to local ULIRGs it is not surprising that we did 
not detect th e 3.3 u m PAH feature. We note however 
that llmanishil (|2006( ) used a 0.9" wide slit for the L-band 
spectroscopy, and the IRS Short-Low slit is 3.7" wide. 
The nine ULIRGs that we use here are 2-4" across so 
there is a concern that the L-band spectroscopy is suffer- 
ing more slit loss than the Spitzer observations. There- 
fore, if the spatial extent of the PAH emitting region is 
wider than the 0.9" slit width, then the total Lq ^/Lj,^ 
ratios will be lower than the derived values, and may 
t hen disagree with the ratios measured in cB58. 

iMoorwoodI (|1986( l measured 3.3 /im line fluxes for a 
number of starburst and Seyfert galaxies and compared 
them to the 8.6 /im fluxes available for seven of them. 
These gal axies have lower lumi nosities than the ULIRG 
sample of llmanishi et al.l (j2007l ). The average luminosity 
ratio, Ls.e/Ls.s = 2.1, is higher than our 3tT lower limit 
of L8.6/-^^3.3(cB58) > 1.9 (as are 4 of the 7 individual 
ratios). Thus, our limit is compatible with many of these 
lower lumin osity galaxies a s well a s the ULIRGs. 

Recently, iMagneUi etal] (|2008f) flnd evide nee for 3.3 
/im PAH emission in 5 LIRGs at 0.638 < z < 0.839. In 
these flve LIRGs they give a range of line-to-L/i? ratios 
L3.3/L1R < 1.29 - 3.5 X 10-3. The 3a 3.3 /im limit 
for cB58 corresponds to L^^^/Lin < 0.003, so we would 
expect marginal (2-3cr) detections if the ratio were similar 
in cB58. 

We also do not detect the Paa line at Xrest — 1.87/im. 
If we assume Case B recombination and scale from the 
Ha flux with a reddening correction using the Calzetti 
re ddening law a nd th e gas reddening correction derived 
m iBaker et all (|2001h . we find that the expected Paa 
flux, fpaa = 2 X IQ-'^^ ergs s"'^ cm"^. This corresponds 
to about a la signal if all of the flux is contained within 
one pixel. Thus, it is not surprising that we do not detect 
the Paa. 

4. DISCUSSION 

4.1. Star Formation Rate and Reddening 

There are now three independent indicators of the star 
formation rate of cB58: rest-frame UV continuum, HI 
recombination lines, and infrared luminosity. When con- 
verting these measurements to star formation rates, we 
assume a Sal p eter I MF with a 0.1-100 M© mass range. 
iPettini et al.l (|2000D . adopting an LMC reddening law 
(|FitzpatricklU986f) . estimate a UV-derived SFR{UV) = 
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93 ± 11 Mq yr^i (with a Salpter IMF between 0.1-100 
Mq). Continuous star formation for 100 Myr was as- 
sumed, but if a more realistic younger age of 9 Myr 
is assumed, the rate must be adjusted higher by 57% 
(|Kenniciittlll998D to SFR{UV) = 146±18 M© yr"!. Our 
estimate based on the UV/optical/IRAC SED assuming 
a Calzetti reddening law is SFR{UV) - 98^;^^ Mq yr"!. 
The Ha line flux, corrected using a Calzetti reddening 
law and a color excess of the g as of E{B — V)gas = 0.06 
derived from the iTeplitz et al.l |20001 ) I3almcr decrement 
of Ha:H/3 = 3.09, giv es SFR(Ha) = 2 5 ± 1 Mq yr"! us- 
ing the conversion of lKennicutlj (|1998f ). Finally, we con- 
vert our infrared luminosity to star formation rate using 
iKennicutli (fl998l ). yielding SFR{IR) = 19 Mq yr"! and 
a maximum of SFR{IR) < 41 Mq yr-i (assuming the 
best-fit DH02 model). The derived SFR derivations are 
summarized in Table [21 

The Ha and IR-derived SFRs agree with each other, 
while the UV-derived SFR is higher by a factor of 3 — 5. 
Of course, the IR SFR assumes optically thick dust, such 
that the vast majority of the glaxy's bolometric lumi- 
nosity is emitted in the IR. For a UV-luminous galaxy 
such as cB58, we should consider the uncorrected UV 
star-formation rate (ie. using observed UV luminosity) 
in addition to the IR-derived SFR. The uncorrected UV 
SFR, hsted in Table [2] is S F R(UV„,n.r.^rr) = 17M q j/r"! 
when using the conversion of lKennicutti (|1998l ). Thus, 
even the addition of this component gives us a total SFR 
a factor of 2-4 times lower than the UV-derived SFRs. 
Furthermore, cB58 is not significantly detected in a deep 
(50ks) Chandra image, also suggesting that the SFR 
must be lower than the UV-derived SFR (O. Almaini 
2008, private communication). It has been suggested 
previously that the discrepancy between the IR- and UV- 
derived SFRs is due to an overestim ate of the redden- 
ing based on t he UV continuum slope (|Baker et al.ll200ll : 
ISawickil [2001h . However, these previous L/_r estimates 
were highly uncertain and could have been underesti- 
mated if the dust was significantly warmer than in local 
starbursts. Our IR data at shorter wavelengths confirms 
that the dust tempera ture can not be m uch more than 
the 33 K assumed bv iBaker et al.l ()2001[ ). In fact, our 
well constrained maximum L/jj is still significantly lower 
than the UV SFR suggests it should be. Therefore, the 
Calzetti reddening law is likely implying too high a UV 
ex tinction for cB58. T his is consistent with the analysis 
of iReddv et"al] (|2006D . who find that very young LBGs 
(<100 Myr) typically have UV slopes that overpredict 
their IR lumi n ositie s when using the relation derived by 
iMeurer et al.l ()1999[ ). suggesting that a different extinc- 
tion curve may be appropriate for young LBGs. 

cB58 exhibits a damp ed Lya profile, indicative o f 
high HI column densities (jSchaerer fc Verhammd [20081) . 
and opaque low ionization metal absorption lines sug- 
gesting a covering fraction of neutral gas (and presum- 
ably dust) near unity (Pcttini et al. 2002; Shaplcy et alj 
l2003f) . Therefore, a reddening law for dust distributed in 
a uniform foreground screen may best describe the com- 
bined UV/IR properties of cB58. The Calzetti reddening 
law is an empirical relation from many starburst galaxies 
with varying dust geometries, and is strongly weighted by 
the relatively unobscurred star-forming regions in these 
systems. Thus, the Calzetti law may be more "grey" 



than is appropriate for galaxies with such large covering 
fractions of gas and dust, such as cB58. When using 
line-of-sight ext inction curves that are steeper in the UV 
fSMC fc LMC. lPrevot et all [1981 [FiFzpatrickl[l986l . re- 
spectively) to fit the stellar SED, we derive smaller star 
formation rates that are closer to the Ha, IR, and X- 
ray SFR derivations. In fact, if we take out the 2175 A 
feature from the LMC extinction curve, we get a much 
better fit to the photometry than when using the Calzetti 
law (reduced = 1-2 and 1.7, respectively), as it repro- 
duces the steep UV spectral slope (see inset in Figure 
[2]) . The SFR derived when using LMC extinction (with- 
out the 2175 A feature) is 2.6 times lower than with a 
Calzetti law {SFRlmc{cB58) = 38 MQyr^^, see Table 
[2| for a summary) and more in line with the other SFR 
determinations. The better fit and realistic SFR when 
using the LMC curve, suggests that a steeper curve than 
the Calzetti law may be more appropriate for cB58. 

This physical interpretation of a large covering frac- 
tion of neutral gas and dust in the foreground of cB58 
may also help explain the discrepancies between UV- 
and I R-derived SFRs in other young {tage < 100 Myr) 
LBGs (jReddv et al.ll2006f l. The spectra of the youngest 
LBGs, like cB58, also exhibit higher equivalent width 
low-ionization interstellar absorption lines and redder 
UV spectral slopes than their older counterparts. If the 
higher equivalent width absorption lines in these systems 
can be attributed to higher covering fractions, rather 
than higher velocity dispersions (the velocity structure 
is not resolved in the low-resolution spectra of unlensed 
LBGs), then a uniform foreground sheet dust geome- 
try (and associated steeper extinction curves) may ex- 
plain the relatively red UV spectral slopes of these young 
LBGs. 

It is possible that the high frequency of Type II super- 
novae (SNe) in a young, intense, and compact starburst 
such as cB58 can significantly alter the composition of 
the dust from what is measured in more quiescent and 
older star-forming galaxies in the local universe (where 
extinction curves have been measured) . The extinction 
curves for Type II SNe dust ([Maiolino et al.ir2004 ). dc- 
rived from dust compos itions expected in Type II SNe 
(jTodini fc Ferraral 1200 1[ ) have been invoked to explain 
the spectrum of a ve ry high redshift BAL QSO (z=6.2, 
iMaiolino et all f200?) a nd the UV SEP of a z = 6.29 
GRB afterglow (Stratt a et al.l [20071 ). These systems are 
very young so their dust is therefore assurn e d to orig- 
inate in Type II SNe. The IMaiolino etld] (|2004[ ) ex- 
tinction curve is flat in the near-UV, but very steep at 
A < 1800 A. Using such an extinction curve, a red UV 
spectral slope between 1300 A < A < 1800 A implies less 
total dust absorption over the entire optical-UV spec- 
trum than would be inferred with an SMC-like extinction 
curve or the Calzetti starburst reddening law. Thus, an 
extinction curve of this nature would help explain the 
lower-than-expected Lu^/Luv because there is less total 
dust absorption, and would help explain the lower-than- 
expected Lho/Luv because the derived intrinsic far-UV 
flux (and thus Lho) wo uld be lower. Howeve r, our fit to 
the steflar SED with the lMaiolino et all (|2004f ) extinction 
curve, though nicely fitting the red UV spectral slope, is 
a poor fit (reduced = 3.6) relative to the LMC and 
Calzetti fits as it does not fit the near-IR (rest-frame 
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optical) photometry. 

A different IMF may also help explain the discrepancy 
in SFR indicators. We have assumed a Salpeter IMF 
that extends up to 100 Mq. If the IMF were deficient in 
the most massive O and B stars, the far-UV flux would 
be much lower than previously predicted, resulting in 
lower-than-expected Lhu and Lip. However, the rest- 
frame UV spectrum of cB58 is fit very well with stellar 
population models with a Salpeter IMF and an up per 
mass limit Mup = 100 M0 (jPettini et al.ll2000ll2002l ). In 
particular, the observed P-Cygni profile of C IV requires 
the p resence of stars with M > 50 Mq (jPettini et al.l 
120001) . Therefore, it is unlikely that a markedly different 
IMF is the cause of the discrepancy between UV- and 
IR/Ha-derived SFRs. 

4.2. Infrared Comparison with Other Starbursts 

4.2.1. Mid-IR Spectra Comparison 

In the upper panel of Figure IH we plot the infrared 
spectrum of the co re of NGC 7714 (Lip = 5.6 x lO^^Lg, 
iBrandl et al.ll200"^ and a composite spectrum of 13 local 
starbursts {<log{Lip)>= 10.7 Lg, [Bran dl et al . 200d), 
scaled to match the flux of the 6.2 /xm feature. The IR 
luminosity of cB58 is 2-4 t imes higher than th at of NGC 
7714 and the starbursts of iBrandl et~aLl (j2006f ). The rel- 
ative fluxes of the PAH features of cB58 agree very well 
with those seen in NGC 7714. The iy. 7/^6.2 ratio is 
a little higher in the starburst composite spectrum, but 
this makes little difference to the total 6-9 /xm flux. The 
PAH equivalent widths of cB58 are as large as in the lo- 
cal starbursts. The exact continuum level is difficult to 
discern because of the large number of emission features 
and the limited wavelength coverage, but it is clear that 
the 6-9 /im luminosity is dominated by PAH emission. 

In the lower panel of Figure [HI we compare the mid-IR 
spectrum of cB58 with high redshift, higher luminosity 
starbursts: a composit e spect rum of 12 submm-sclected 
galaxies bv lPope et al.l (j2008[ ) and a lensed, highly mag- 
nified submm source at z = 2. 516 with L ~ 8 x 10^^ 
L© (SMM J163554.2-h661225, iRigbv et al.ll2008l) . AU 
templates are normalized to match the 6.2 fj,m emission 
strength. The submm-selected galaxies are 1-2 orders 
of magnitude more luminous than cB58, and are dustier 
than UV-selected LBGs. The mid-IR spectra look very 
similar both in their line strengths and the shape of the 
continuum, but few small differences can be seen. First, 
as seen in comparing to the Brandl et al. (2006) compos- 
ite of local starbursts, the £7.7/^^ 6.2 is lower i n cB58 than 
in the other spectra. Second, the lPope et al.l (j2008) com- 
posite has a higher continuum to line ratio (especially at 
6.2 /im) than that of cB58, likely due to hot dust emis- 
sion from AGN, whic h may contribute up to 30% of the 
mid-IR fiux in SMGs (iPope et al.ll2"0Cil . 

4.2.2. Mid- to Far-IR Comparison 

The PAH luminosities of pure starbursts have been 
shown to correlate with t otal infrared luminosity (and 



thus star formation rate) (iForster Schreiber et al 



Peeters et all I2Q04I: IBrandl et al.ll2006l: iPope et al 



2004 



2008 



Rigbv et al.l 120081 ). though with considerable scatter. 



This correlation has not yet been tested for high redshift 
galaxies with L < 10^^ Lq such as cB58. We c ompare the 
LpAH and Lip of cB58 to the relations fit bv lPope et al.l 



(j2008f) to local starbursts (jBrandl et al.ll2006f) and high-z 
submm-selected galaxies. First, we re-measure the P AH 
luminosities in the same manner as iPope et al.l (j2008D by 
fixing the continuum levels at either side of the strong 
PAH features (as opposed to simultaneously fitting all 
PAH features as well as the continuum). The assumed 
continuum is shown as the dot-dashed line in Figure [5] 
and the corresponding luminosities are listed in Table |3l 
In Figure [3 which shows the comparison between 6.2 /im 
luminosity and hip, we plot the cB58 position relative 
to the local starbursts and high-z SMGs. The PAH lu- 
minosities lie on the high side of the measured relation 
such that the Lpah is high given the known Lip, but 
the offset from the relation is within the scatter. The 7.7 
/im PAH luminosity is also hig her than the fi tted rela- 
tion by the same factor. The Pope et al.l (j2008f ) relations 
for both the 6.2 and 7.7 /im features both suggest an 
Lip ~ 3 X 10^^ Lq, about a factor of two higher than 
the best estimate of L7i^(cB58). 

Another important relation is the correl ation of the 
rest-frame 8 /im and total IR luminosities (jRigbv et al.l 
120081 iBavouzet et al.l [20081: ISaiina et al.l I2008D . This is 
particularly useful for broadband only studies at high 
redshift (eg. Spitzer 24 /zm) where individual PAH lu- 
minosities can not easily be measured. From the IRS 
spectrum, we measure the rest-frame 8 fim luminosity 
(the luminosity measured through the IRAC 8 /im band 
at ^ = 0) of cB58 to be Lf^(cB bi) = 2.57 x 10^° L©. 
Using the fitted relation of Ba vouzet et al.l ([2008), the 
L8(c-B58) predicts an Lip{cB58) = 1.7 x 10" Lq, in 
good agreement with our Lip estimate. 

Finally, we compare the /24/ .f85Q flux ratio of cB58 to 
that of SMGs (jPope et al.ll20()6f) at similar redshifts. The 
cB58 ratio, /24//850 = 0.057, is higher than six of the 
seven SMGs at similar redshifts (2.3 < z < 3.1), which 
range from 0.0031 < 124,/ f 850 (SMGs) < 0.062. The one 
SMG with a similar ratio has an X-ray hardness ratio in- 
dicative of an obscurred AGN th at is likely c o ntirib uting 
significantly to the 24 /im flux. iPope et all (j2006D sug- 
gest that additional extinction of the rest-frame mid-IR 
flux is required to reproduce the low /24//850 flux ra- 
tios observed in SMGs. Therefore, the higher /24//850 of 
cB58 may simply be due to the lack of extinction in the 
rest-frame mid-IR relative to SMGs. Additional extinc- 
tion of SMGs in the mid-IR is expected, as these systems 
contain far more dust than LBGs and often e xhibit severe 
extincion in the rest-frame UV/optical (eg. iSmail et al.l 
Il999i) . Alternatively, the discrepancy can be explained 
by more cold dust emission in SMGs relative to cB58. 
This is perhaps not surprising as SMGs are selected on 
their cold dust emission. 

5. CONCLUSIONS 

We have obtained Spitzer photometry and mid-IR 
spectroscopy of the lensed LBG, MS1512-cB58. This 
is the (intrinsically) least luminous galaxy to be stud- 
ied in detail in the infrared at 2 > 1. Although cB58 
is a "typical" LBG in terms of its UV and IR lumi- 
nosities, the lack of a Balmer break in the SED indi- 
cates that it is much younger {tage ~ lOMyr) than a 
large majority of L BGs. As seen in previous studies 
(jShaplev et al.l l2005'). stellar population fits to the com- 
bined optical/near- IR/IRAC photometry give similar pa- 
rameters to fits to the optical/near-IR photometry alone. 
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but with smaller error bars. The mid-IR photometry do 
not significantly constrain the mass of an older stellar 
population, as the rest-frame near IR flux from the cur- 
rent burst {Mburst ~ 10® M0) can dwarf the flux from 
an older stellar population with twice the mass. 

The far-IR photometry is reasonably well fit by local 
starburst templates and gives an Lm = 1 — 2 x 10^^ L0. 
In addition, the PAH luminosities and rest-frame 8 /im 
luminosities agree with the LpAfj-io-Ljn and L^-to-Ljn 
relations measured in local starbursts of comparable lu- 
minosity and high redshift, higher luminosity starbursts. 

The inferred SFR from the infrared luminosity is con- 
sistent with the SFR derived from the extinction cor- 
rected Ha luminosity, but the UV-derived (extinction 
corrected) SFR is larger by a factor of 3 -5. This phe- 
nomenon has been noted bv iReddv et al.l ([200 6). where 
LBGs with starburst ages less than 100 Myr consistently 
have lower Lm than the UV spectral slope and lumi- 
nosity suggest. With cB58, we can be confident that 
this discrepancy is not due to a poor determination of 
the infrared luminosity. This suggests that the Calzetti 
obscuration law for starbursts may not be valid in very 
young, < 100 Myr old, high redshift starbursts. 

The high HI column density and large covering fraction 
of neutral gas suggests a dust geometry approximated by 
a uniform foreground sheet. This geometry would result 
in a steeper reddening law like the LMC or SMC curves 
and would explain the overstimate of the UV obscuration 
using the Calzetti reddening law. In addition, because 
many young {tage < 100 Myr) LBGs also exhibit spectral 
properites indicative of large covering fractions of neu- 
tral gas (eg. high equivalent width low ionization metal 
absorption lines), a steeper reddening law may also be 
appropriate for explaining the relatively red UV spectral 
slopes observed in these systems as well. 

It is of course impossible to apply the relations mea- 
sured in cB58 to the LBG population as a whole. For 
example, there is significant scatter (beyond measure- 
ment errors) in the predicted opacities versus UV spec- 
tral slope in the starburst sample used to derive the 
Calzetti law (Calzetti et al. 1994). Therefore, cB58 may 
simply lie to one end of the natural dispersion already ob- 
served in low redshift starbursts. However, the fact that 
the phenomenon of the UV overpredicting Ljr has been 
observed in a large sample of young LBGs (|Reddv et al.l 
[2OO6.) suggests that cB58 may be typical of these young 
sytems, and that a systematic bias in UV-derived SFR 
derivations may exist for the youngest LBGs. SED fits to 
rest-frame UV and optical photometry of LBGs suggest 
that at least ^ 30% hav e starburst ages less than 100 Myr 
(jShaplev et all 2001; Pa povich et al.ll2001h and, at least 
for LBGs with Ljjv ^ L* , there is no correlation between 
age and observed UV luminosity (|Shaplev et al.„ .200l]). 
Therefore, if the SFR in all of these young LBGs is over- 
estimated by a factor of ~ 4 — 5, then the global star for- 
mation rate density of all LBGs will be overestimated by 
a factor of 2. S pitzer obs ervations of recently discov- 
ered lensed L BGs (iSma fl et al. 2007; Coppin et al.ll2007l : 
lAllam et al.l l2007f ) will allow us to better determine if 
these discrepancies persist in the rest of the young LBG 
population. 
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TABLE 1 

Spitzer Photometry of MS1512-CB58 



Band 


Flux Density 


Error 








O.O 


oy.u 




4.5 


77.5 


5.5 


5.8 


77.0 


4.6 


8.0 


53.7 


3.8 


16 


87 


10 


24 


240 


40. 


70 


1.7 mjy 


1.0 mJy 


160^ 


<24 mJy 




3a confusion limit fromlDole et al.l l|2004l) 
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TABLE 2 
Star Formation Rate Estimates 



SFR Indicator SFR jv 
Mq yr-i 

UV 

SED Fit (This work) 

Calzetti Law QSlg^ 1.7 

LMC 38lig 1.2 

Pettini et aL 2000*' 146 ± 18=^ 

Uncorrected -Lieoo 17 it 1'^ 

IR (This work) 19-41= 

Hg'^ 25 ± V 

Corrected for 30x magnification'' Corrected to use Ljjy-to-SFR relation for age of 9 Myr, rather than 100 Myr'^ Errors are statistical 
and do not reflect uncertainties in conversion to SFRs'^ Extinction corrected assuming E(B-V)gas = 0.12 UTepIitz et al.ll2000l : iBaker et al.l 
[20011) 



TABLE 3 
Emission Line Fluxes of cB58 



Wavelength 
[/im] 




Flux 

[lO^-"^^ erg s^^ cm~^] 


Luminosity ^ 
[10** erg s-i] 


Rest Equivalent Width 
[^tm] 


1.87 (Pact) 




<0.60 


<0.37 




3.3 "° 




<0.84 


<0.52 




6.2 




3.3+0.6 


2.0 (1.1)=^ 


5.4 


7.7 (7.42+7.60+7.85) 




10.0+3.0 


5.8 (2.7)== 


12 


8.6 




1.6+0.8 


1.0 


1.6 


Not corrected for ma; 
as a dot-dashed line in 


gnification'' 3cr 
Figure |5] 


limif^ Values in parentheses 


are the line luminosities derived when 


assuming the continuum shown 




Fig. 1.— HST/IRAC/MIPS postage stamps of MS1512-cB58. From top left: HST NICMOS F160W (H-band) followed by the four IRAC 
bands, IRS 16 /^m, MIPS 24 and 70 /xm. cB58 is marked with a circle of 3" radius, and the cluster cD galaxy is marked with a white cross. 
The cluster cD and cB58 are ~ 5.3" apart. The slit positions and position angles are shown in the images at the appropriate wavelengths 
(Short-Low in IRAC 5.8 & 8.0 /jm, and Long-Low in IRS 16 and MIPS 24). The Long-Low slit is offset 2.35" from cB58 in the direction 
opposite from the cD galaxy (i.e. 27.5° W of N). The postage stamps are 34"x41" and North is up. East is left. 
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Fig. 2. — The SED fit to the optical/near-IR/mid-IR photometry of cB58. The sohd thick blue line is the single component fit using a 
Calzetti obscuration law. The thin solid line is the sum of two components: a continuous star formation component (dashed) and an old (2 
Gyr) single burst component (dotted). The combined two component SED is similar to the single component SED and is therefore hidden 
under the thick blue solid line. Also plotted (green dot-dashed line) is the single component fit with an LMC extinction curve (without the 
2175 A feature) that better produces the red UV spectral slope (shown in the inset). 
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Fig. 3. — The infrared SED of cB58. The MIPS 16, 24 and 70 /im detections are plotted along with the conservati ve 70 and 160 Atrn 
limits. The 850 fivci detection and limit ( van dor Werf et al. 2001; Sawicki 2001, respectively) and the 1.2 mm detection IIBaker et al.ll200]I) 
are also plotted. The best-fit SED from lCharv fc Elbaj ll20dlfl is plotted (solid line), giving a total Lm = 1.1 X 10^^ L(n . The d ashed curve 
is the CEOl SED when picking a local template based on the L24, a lone, the dash-dotted curve is the best fit SED of lDale He He lou (200!^), 
and the dotted curve is the IR SED assumed bv lBaker et al.l I I2001I) . The first-order Long- Low spectrum is shown in blue for reference. 
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Fig. 4. — The IRS spectrum of cB58. Short-Low2, Short-Lowl, Long-Low2, and Long-Lowl are plotted in purple, blue, green, and 
red solid lines, respectively (errors are plotted as dashed lines). The wavelengths of Paa and the primary PAH emission features are 
labeled. The original Long-Lowl spectrum before corrections for cD galaxy contamination and slit loss is plotted (thin solid grey line) to 
demonstrate how little these corrections affect the spectrum. 
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Fig. 5. — Smoothed Long-Low 1st order spectrum of cB58 (black histogram) and the best-fit spectrum (red). The best-fit indi vidual PAH 
profile s (blue solid lines) and continuum (green dashed line) are also plotted. The assumed continuum when comparing to the lPope et al.l 
II2008I ) measurements is shown as a magenta dot-dashed line. 
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Fig. 6. — Comparison of the mid-IR spectra of cB58 (smoothed, black histogram) to other starbursts. All comparison spectra are 
normalized to have the same p eak flux at 6.2 iim. In the upper panel we compare to compa rably bright l ocal starbursts: the spectrum 
of NGC 7714 (blue dashed line, IBrandl et al.l l2004l ) and the average starburst spectrum from Bran dl et al.l l|2006. red dotted line). In the 
lower panel we compare t o high re dsliift starbursts: a composite spectrum of SMGs more than an order of magnitude more luminous than 
cB58 (purple dashed line. lPope et al.ll2008l) . and a lensed starburst with 4 — 6 times the Lju of cB58 (green dotted line. lRigbv et al.ii2008i ). 
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Fig. 7. — The measure d correlatio n between 6.2 /im and total infrared luminosities for local starbursts'Brand l et al.l ( open square s . 120061) 
and high redshift SMGs lPope et aP (open diamonds, 2008). This figure is adapted from the top panel of Figure 12 in lPope et al.l 1)20081 ) . 
The best-fit relation to both populations is plotted as a dashed line. cB58 (large open circle) lies above the relation such that it has more 
PAH luminosity relative to its infrared luminosity, but it lies within the scatter of the other starbursts. 



